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Non-equilibrium bosonization technique facilitates the solution of a number of important many- 
body problems out of equilibrium, including the Fermi-edge singularity, the tunneling spectroscopy 
and full counting statistics of interacting fermions forming a Luttinger liquid. We generalize the 
method to non-equilibrium hard-core bosons (Tonks-Girardeau gas) and establish interrelations 
between all these problems. The results can be expressed in terms of Fredholm determinants of 
Toeplitz type. We analyze the long time asymptotics of such determinants, using Szego and Fisher- 
Hartwig theorems. Our analysis yields dephasing rates as well as power-law scaling behavior, with 
exponents depending not only on the interaction strength but also on the non-equilibrium state of 
the system. 
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PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 



I. INTRODUCTION 

There is a number of quantum many-body problems 
that are of paramount importance for condensed mat- 
ter physics and, at the same time, possess an exact solu- 
tion. These are the Anderson orthogonality catastrophe^, 
Fermi edge sing ularity 2 - (FES), Luttinger liquid^ (LL) 
zero-bias anomaly^, and Kondo problems^. It has been 
realized long ago that these problems are, in fact, inter- 
connected, both vis-a-vis the underlying physics and the 
mathematics involved. Such connections have been used, 
e.g., for the representation of the dynamics of the Kondo 
problem as an infinite sequence of Fermi-edge-singularity 
events^. 

In our recent works^ we have addressed some of these 
problems away from equilibrium. To achieve this goal, we 
have developed a non-equilibrium bosonization technique 
generalizing the conventional bosonization^— onto prob- 
lems with non-equilibrium distribution functions. We 
have shown that the relevant correlation functions can be 
expressed through Fredholm determinants of a "count- 
ing" operator. The information on the specific type of 
the problem, as well as on different aspects of the in- 
teraction, is encoded in the time-dependent scattering 
phase of the counting operator. The results of Refs. 0H 
have demonstrated that those classic many-body prob- 
lems are more closely connected than has been previ- 
ously understood, extending the interrelations into the 
non-equilibrium regime. 

The purpose of the present paper is twofold. First, 
we present the relation between the LL tunneling spec- 
troscopy, the full counting statistics, and the Fermi-edge 
singularity problems in a systematic way. We also in- 
clude in this scheme the generalization of the LL prob- 
lem to a non-equilibrium gas of strongly repulsive bosons 
(Tonks-Girardeau gas). These results are summarized in 
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TABLE I: Non-equilibrium correlation functions of many- 
body problems: Fermi edge singularity (Gp ES ), Green func- 
tions of right- and left-moving fermions in a LL (Gp R and 
Gp L ), full counting statistics of a LL (xr(A)), Green func- 
tion of the Tonks-Girardeau gas (Gg). All these correlation 
functions can be cast in the form G > = (e- iO - (T) e iO + (0) ), 
where O = ^2 v=a L c v <j) v and 4> v are free bosonic fields. The 
coefficients c v are shown in the second and third columns of 
the table. Evaluation of these correlation functions yields 
the results in the form of Fredholm- Toeplitz determinants 
A_r[i5_r] Al [Sl]- The corresponding phases <5r, l are presented 
in the last two columns for the case where the LL is adiabat- 
ically connected to the reservoirs. 



Table |IJ Second, we apply the general theory of asymp- 
totic behavior of Toeplitz determinants to the present 
problems. This allows us to extract important informa- 
tion on the many-body physics of these models, including 
non-equilibrium dephasing rates and the modification of 
power-law exponents due to non-equilibrium conditions. 

The structure of the present paper is the following: In 
Sec. Ill Al we summarize the results of our earlier works^ 
for the FES problem, as well as the tunneling spec- 
troscopy and the counting statistics problems in the LL, 
all underlined by interacting fermions. Section fllBl is de- 
voted to the tunneling spectroscopy problem for a non- 
equilibrium Tonks-Girardeau gas. In Sec. IIIII wc study 
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the long-time behavior of emerging Fredholom determi- 
nants, employing Szego and Fisher-Hartwig theorems. 
For completeness, properties of the Tonks-Girardcau gas 
at equilibrium are summarized in Appendix^ while rele- 
vant mathematical theorems from the theory of Fredholm 
determinants are presented in Appendix [Bj 



II. MANY-BODY PROBLEMS AS FREDHOLM 
DETERMINANTS 

A. Fermionic problems 



Free fermions 



We first discuss non-interacting electrons. Evidently, 
the Green functions of free fermions Gq , Gq (we use the 
standard notations of Kcldysh formalism^ can easily be 
calculated within a fermionic description. However, it is 
insightful to recalculate it employing the non-equilibrium 
bosonization- 7 -. The result reads 
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2ttv t =p i/A 



A„[M*)] 



(1) 



Here 77 is the chirality index (77 = L for left and 77 = R 
for right movers), the phase S v (t) = Aw T (i, 0), where 



Writ, t) = 6(t -t)-6(t-t-r) 



(2) 



is a "window function", and A = 2ir. Further, A^[^(/j)] 
is the Fredholm determinant of a counting operator 



A„ [<$„(*)] = det[l + (e-^ - l)n n ] 



(3) 



and we denote by A,, the determinant normalized to its 
value for zero-T equilibrium distribution. As we see, in 
the bosonic description, a free fermion is represented by a 
2ir phase soliton entering the counting operator. The lat- 
ter consists of a fermionic distribution function n v (e) , and 
a time-dependent scattering phase 5 n (t), with e and t to 
be understood as canonically conjugate variables. Since 
the matrices 5 and h can not be diagonalized simulta- 
neously, the task of calculating A 7) [J 7) (f)] is non-trivial. 
In most cases the analysis of such determinants can be 
done only numerically. However, the long time asymp- 
totic properties can be studied analytically, by means 
of Szego and Fisher-Hartwig theorems, as discussed in 

sec.nni 

When read from right to left, Eq. (TTJ) can be viewed as a 
remarkable identity for Fredholm determinants. It yields 
an exact value of such a determinant for an arbitrary 
distribution function n n (e), and for the phase function 
5(t) given by a window function with amplitude 2n. It 
is worth mentioning that the determinat should be con- 
sidered as analytically continued from the region of small 
phases 6. 



2. Fermi edge singularity 

The FES problem describes the scattering of conduc- 
tion electrons off a localized hole, which is left behind 
by an electron excited into the conduction band. Histor- 
ically, the FES problem was first solved by exact sum- 
mation of an infinite diagrammatic series^. While in the 
FES problem there is no interaction between electrons in 
the conducting band, it has many features characteris- 
tic of genuine many-body physics. Despite the fact that 
conventional experimental realizations of FES are three- 
dimensional, the problem can be reduced (due to the local 
and isotropic character of the interaction with the core 
hole) to that of one-dimensional chiral fermions. For this 
reason, bosonization technique can be effectively applied, 
leading to an alternative and very elegant solution^. 

One can consider the FES out of equilibrium^, with 
an arbitrary electron distribution function n(e). This 
problem can be solved within the framework of non- 
equilibrium bosonization^, with the following results for 
the emission/absorption rates: 



G FEs( T ) 



iAA t (2tt- 28q) 
' 2ttv(1 ± iAT^-^M 2 



(4) 



Here So is the s-wave electronic phase shift due to the 
scattering of conduction electrons off the core hole. 



3. Luttinger liquid: Tunneling spectroscopy 

The tunneling spectroscopy technique allows one to ex- 
plore experimentally Keldysh Green functions of an inter- 
acting system that carry information about both tunnel- 
ing density of states and energy distribution. Recent ex- 
periments on carbon nanotubes and quantum Hall edges 
have proved the efficiency of this technique in the con- 
text of ID system o 17 ' 18 . The technological and experi- 
mental advances motivate the theoretical interest in the 
tunneling spectroscopy of strongly correlated ID struc- 
tures away from equilibrium£ii2r— . 

In the case of a LL made of ID interacting fermions, the 
Keldysh Green function may be evaluated theoretically 
via the non-equilibrium bosonization technique. Assum- 
ing that a long LL conductor is coupled to two reservoirs 
(modeled as non-interacting ID wires^— ) with distribu- 
tion functions nn(e) and respectively, one obtains 7 - 

3 _^ iA A R [S R (t)]A L [5 L (t)} 
Gr{t) ~ T 2^u (l±zAr)i+7 ' (5) 

where u = v/K is the sound velocity, 

1 =(1-K) 2 /2K, (6) 

and 

K = (1 + g/TTv)~ 1/2 (7) 
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is the standard LL parameter in the interacting region. 
The phase S v (t) is found to be a superposition of rectan- 
gular pulses, 

oo 

^7 ft) = ^ ^v,nW T (t, t n ) , (8) 
n=0 

where 

t n = (n + 1/2 - l/2K)L/u (9) 

and 

S v ,2m = nt- n rZr%(l + r,K)/VK , 

^, 2m+ i = -Trt-^+V™ (1 - v K)/Vk . (10) 

Here r,, , t v are reflection and transmission coefficients of 
plasmons at the left (77 = L = — 1) and right (rj = R = 
+1) boundaries. For r <C L/it the coherence of plasmon 
scattering may be neglected and the result assumes the 
form of a product 

00 

A v [S v {t)} ~ Y[ A VT (S v . n ). (11) 

n=0 

The plasmon scattering causes a fractionalization (c/. 
Refs. l27lj29l - [33 ) of the phase soliton, splitting it into an 
infinite series of pulses. As a result, the Frcdholm de- 
terminant of a counting operator takes the form of an 
infinite product of determinants, each calculated for a 
rectangular pulse with a corresponding scattering phase 

Sr,,n(t) = S, hn W T (t,0). 

In the case of smooth boundaries with the leads we 
have r v ~ 0, which implies that only the first (n = 0) 
pulse survives in both Sr and Sl- The corresponding 
amplitudes are depicted in Table |U 

The above results correspond to the case where the 
tunneling spectroscopy point is located inside the inter- 
acting part of the wire. The case of tunneling into one 
of non-interacting leads is studied in the same way; the 
results can be found in Ref. 0- 

So far we have discussed Green functions at coincid- 
ing spatial points having in mind tunneling spectroscopy 
experiments. Green functions at different spatial coor- 
dinates are also of interest, in particular, in the con- 
text of Aharonov-Bohm interferometryi 3 -^. (A similar 
problem in the context of chiral edge states has been 
considered in Refs. I3514371 .) Equation ([S]) expressing the 
non-equilibrium LL Green function in terms of Frcdholm 
determinants can be generalized for this case as well^. 

4- Luttinger liquid: Counting statistics 

In the problem of counting statistics of non-equilibrium 
LL one is interested in the generating function xW — 
^^° = _ 00 p T (n)e mA , where p T {n) is the probability for n 
electrons to pass through a given cross-section during 



the time interval r. On the experimental side, the sec- 
ond moment (shot noise) has been measured in correlated 
ID system s 38 i 39 ; an experimental analysis of higher mo- 
ments and of the full counting statistics (FCS) remains 
a challenging issue. 

For the non-interacting case the FCS generating func- 
tion xW nas been calculated in Ref. [40 by means of 
Landaucr scattering-state approach. For an ideal quan- 
tum wire (with no scattering inside the wire) with dis- 
tributions n v (e), the generating function of FCS is given 
by 

x(X)=A R [6 R (t)]A L [S L (t)}, (12) 

with the phases S^ft) = Xrjw T (t,0). For an interacting 
wire (LL) the result^ retains the form of a product of 
Frcdholm determinants, Eq. (fl"2|) . but with the scattering 
phase <5^(i) turning into the following sequence of pulses: 

00 

S v {t) = '^2S ri!n w T (t,t n ) , (13) 

71=0 

with partial phase shifts 

S v ,2n = r)Xt-r,"/Kr™r" v = r)\e* 2n , (14) 
S v ,2n+i = vX-t-n^ r n +lr -ri = V^^n+i ■ (15) 

The time moment t n = (n+ 1/2— 1/2 K)L/u corresponds 
to the beginning of the n-th pulse. 

The result is non-trivial when the measurement inter- 
val r is small compared to the length of the LL conductor, 
t <C L/u. (In the opposite limit the pulses overlap and 
one recovers the results of the FCS of non-interacting 
fermions^.) As for the tunneling spectroscopy prob- 
lem, the determinants then split into a product of de- 
terminants for individual pulses. The FCS of the LL 
then becomes a superposition of FCS of non-interacting 
electrons with fractional charges e* n . For the case of 
smooth boundaries we obtain only one fractional charge, 
e* q = \f~K (Table Q}. I n the opposite limit of sharp 
boundaries, we obtain a sequence of fractional charges of 
the form e\ n = 2K{1 - K) n /(1 + K) n+1 . 

These results describe current fluctuations in the in- 
teracting part of the wire. The FCS measured in the 
non- interacting part (keeping the assumption r <C L/u) 
has a similar structure, but the values of the fractional 
charges are different, see Ref. H 

B. Luttinger liquid of bosons: Tonks-Grirardeau 

gas 

We now consider the problem of a non-equilibrium 
LL formed by bosons with strong repulsion. Interacting 
bosons out of equilibrium attract currently a great deal of 
attention, in particular, in connection with experiments 
on cold atoms^. 

To make the connection between the fermionic and the 
bosonic LL particularly transparent, we will adopt a toy 
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model with the same setup, and with the same behavior 
of the LL interaction constant K{x) as the one assumed 
for the fermionic models considered above. Specifically, 
we will assume that K(x) takes takes the value 1 in the 
reservoirs, and a value K in the central part of the setup. 
For a gas of fermions K = 1 implies the absence of in- 
teraction; now, in the present context, it corresponds to 
hard-core bosons (whose many-body wave function van- 
ishes when the coordinates of two particles coincide); this 
system is known as the Tonks- Girardeau gas^. The 
wave function of this system is related to that of non- 
interacting fermions via the transformation 

ip B {xi, ■ ■ -,xn) = s(x u . . .,x N )ip F (x 1) . . .,X N ) , (16) 

where s(x\, . . . , xn) is a sign factor (±1) counting the 
parity of the number of permutations of coordinates, 

N 

s(xi,...,xn) = JJsgn(xj -Xj). (17) 

i>j 

If the fermionic wave function is real, this simply yields 
ipB = \"4>f\- We will further assume that after this boson- 
to-fermion transformation, the fermions in the reservoirs 
are characterized by distribution functions ri^e), as in 
the above fermionic models. While very natural in the 
fermionic language, this requirement is in general quite 
artificial for bosons. We do not know whether it can be 
realized in an experiment and consider this as a theoret- 
ical toy model. A more realistic situation arises in the 
case of partial non-equilibrium, where each of the reser- 
voirs is at equilibrium but the temperatures of the two 
reservoirs are different. 

We will analyze the single-particle bosonic Green func- 
tions 

G>(x,t) = -i(* B (a:, *)*£((), 0)) , (18) 
G<{x,t) = -i(tf B (p,0)* B (x,t)) . (19) 

(20) 

that carry information about spectral properties (density 
of states and distribution functions) of the system. To 
calculate them, we proceed, in analogy with fermionic 
systems, via the non-equilibrium bosonization technique. 
The term "bosonization" here is, perhaps, not optimal 
since the original system is bosonic to begin with. What 
actually happens is a transformation from the original 
bosonic fields to new bosonic fields </>, 0, the latter 
describing density fluctuations in the system. The orig- 
inal field operator is expressed in term of the new fields 
asi£ 

*b0«0 = V/°o + n(x)j J2 e lm ^\e~ W(X) ■ ( 21 ) 

^ mGeven ' 

Here the field 4>(x) is related to the smeared density 
p(x) = po + Tl(x) (where po is the average density) via 



p(x) — —d x cj)(x)/ir. The bosonic fields (4>,0) satisfy the 
commutation relation 

[0(x),^')] = ysgn(x-x') . (22) 

The bosonization prescription (|2ip may be further sim- 
plified by discarding fast oscillating terms (corresponding 
to to 7^ in Eq. (|21|1 ) and neglecting n(a;) in comparison 
to po in the preexponential factor. We thus obtain 

~ v^ e -*<*> . (23) 

The goal of the discussion above was to present the cor- 
relators needed for G B (x, t) in terms of the bosonic fields 
cj>, 6. We next discuss the pertinent action. 

Since the boson-to-fermion transformation (|16|) pre- 
serves the density | - i/ , | 2 i and in view of our assumption 
about the density matrix in the reservoirs, the Keldysh 
action of the non-equilibrium Tonks-Girardeau gas has 
exactly the same form as for the corresponding problem 
of free fermions^, 

S [p,p] =J2 S o[Pv,Pvl > ( 24 ) 
v 

where p and p are the classical and quantum components 
in Keldysh representation, and the action for excitations 
with a given chirality is 

So \p n ,p v ] = -p v n% 1 p n - i In Z v \x v ] ■ (25) 

The information about the non-equilibrium state of the 
problem is encoded in the infinite sum of vacuum loops 

00 .-n+l 

ilaZ ri [x v } = J2^r^S n , v , (26) 

n=2 

representing a partition function of free fermions subject 
to the external quantum field 

X„=nfV (27) 

Here S„ tTI is the n-th order density cumulant of free 
fermions in the given non-equilibrium state . The right 
and left density components are related to the fields <j) 
and 9 via 

Pn(x) = -^d x <f> n , (28) 

where 

4> R = 6 - 4> , 

^L = + cj). (29) 

In order to find the Green functions of the original bosons 
one therefore needs to calculate the functional integral 

G B {x,t) = -ip J V(l)e iS[p ' p] 
xe (-i/V2)[e[o,o]-e[x,t]+§(o,o)+e(^,t)] 
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and similarly for G < . This functional integral is fully 
analogous to the one we had to evaluate while studying 
the Green functions of free fermions (see Ref. for de- 
tails). The only difference is in the source part, i.e., in 
the exponent that contains terms linear in bosonic fields. 
This difference has a rather transparent meaning: cre- 
ation of the fermion ip^ corresponds to the vertex oper- 
ator % Frj e~ % ^ that generates a soliton (a step-like 
plasmon wave with amplitude 2ir) that propagates in the 
direction rj. By contrast, bosons represented by the op- 
erators $b do not have chirality; the corresponding cre- 
ation operator ^ B cx e~ l6 = e^C^R+'W generates both 
right and left moving waves. This results in two Fredholm 
determinants — one corresponding to the left and another 
to the right reservoir — both with scattering phases tt (as 
opposed to a single determinant with a scattering phase 
2tt for free fermions): 

G b( X ' T ) = -ipO^T-x/vM^L^+x/vi^) 
xe -i^f[sgn{T+x/v)+sgn(T-x/v)] 

Here the numerical prefactor has been fixed by compari- 
son with the equilibrium case, see Appendix [A] At zero 
temperature A r (5) ~ (1 + A 2 t 2 )~ 3 ( S / 27T ) ; which repro- 
duces the well-known x -1 / 2 (for r = 0) behavior corre- 
sponding to the fc -1 / 2 momentum distribution of parti- 
cles in the equilibrium T = Tonks-Girardeau gas. 

Equation (f3~Tj) yields the Green function of the non- 
equilibrium Tonks-Girardeau gas with a LL constant 
K = 1. We now assume that the interaction is differ- 
ent in the central part of the system, so that K{x) there 
is equal to K ^ 1. Performing the analysis in analogy 
with the fcrmionic problem, we find 

G>(x 7 t) = -tp A R [6 R (t)}A L [S L (t)}e-^ s ^ , (32) 

where each of the phases S v (t) is given by the arithmetic 
mean of the corresponding phase for fermionic Green 
functions Gfr and Gfl- In general we obtain again an 
infinite sequence of pulses, as for the fermionic problem. 
In the case of smoothly varying K(x) only the first pulse 
survives in each of the phases S^t); its amplitude is equal 
to 7r/ y/K, as shown in Table Q] 

C. Summary 

Let us briefly summarize the results presented above. 
All the problems considered (Fermi edge singularity, tun- 
neling spectroscopy of interacting fermions, their full 
counting statistics, and spectral properties of interacting 
bosons) can be solved by the non-equilibrium bosoniza- 
tion approach, with the results expressed in terms of 
Fredholm determinants of counting operators. In these 
expressions, all differences between the problems are en- 
coded in the values of scattering phases S v (t). These 
scattering phases consist either of one pulse (for the FES 
problem and LL problem with smooth boundaries), or of 



a sequence of well-separated pulses (for a sufficiently long 
LL sample with sharp boundaries). In the latter case the 
determinant can be split into a product of determinants, 
each of which corresponding to a single pulse. Therefore, 
physical properties of a number of many-body problems 
are governed by the behavior of Fredholm determinants 
with a single phase pulse. The analysis of this behavior 
is presented in the next section. 

It is worth mentioning that there is a vast literature 
on the connection of Fredholm determinants to count- 
ing statistics as well as to quantum and classical inte- 
grable models and free-fermion problems; see, in partic- 
ular, Refs. SIMM 

III. ASYMPTOTIC PROPERTIES OF 
FREDHOLM DETERMINANTS 

A. Ultraviolet regularization and reduction to 
Toeplitz form 

We consider a Fredholm determinant of the counting 
operator for scattering phase exhibiting a single pulse of 
an amplitude 5 and duration r, 

5{t) = S x w T (t,0) . (33) 

In this case the Fredholm determinant © is of the 
Toeplitz type. To show this, one defines a projection 
operator P that acts on a function y(t) by restricting it 
to the time interval [0, r]: 

f y(t), for t £ [0,r] 
Py(t) = ! (34) 
I , otherwise . 

For the single-pulse scattering phase (|33p Eq. ^ can 
then be rewritten in the form 

A[6(t)] = det[l + P(e- lS - l)n] . (35) 

Since P 2 = P, we can bring the determinant ([33]) to the 
form 

A[6(t)} = det[l + P{e- lS - l)nP] . (36) 

This determinant still requires an ultraviolet regulariza- 
tion. A possible way to introduce it is to discretize the 
coordinate t by introducing an elementary unit of the 
size At = 7r/A, such that tj = jAt. This corresponds to 
restricting the energy e variable to the range [—A, A]. In 
the formulations of Szego and Hartwig-Fishcr theorems 
that will be applied below the function f(z), that gener- 
ates the matrix is defined in the complex plane of variable 
z, and its values on the unit circle z = e lS parametrized 
by the angle 8 £ [— tt,tt] are important. In Eq. (|36|) this 
function is /(e) = 1 + n(e)(e~ lS — 1). The correspon- 
dence between e and 6 is established by rescaling the en- 
ergy 7re/A = 8. Further, we need to eliminate the jump 
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in n(9) at 9 = ±ir that results from a hard cutoff and 
would generate an additional, unphysical contribution of 
the Fermi-edge type. This is done by introducing a phase 
factor, 



f(e) = [l + n(e)(e-* s - l)]e 



(37) 



that makes /(e) periodic on [—A, A]. Fourier trans- 
forming the periodic function /(e), we obtain f(tj) with 
tj = jV/A. Equation Q36p then reduces to a determinant 
of a large (TV x TV, where TV = tK/tt) but finite matrix 



A N [f] = dot [/fo - t k )} , < j, k < TV - 1 



(38) 



In the above derivation we assumed r > 0; the result for 
t < follows from the property A_ T (J) = A T (— <5). 

The matrix {f(tj - t k )} with < j, fc < TV — 1 is 
of a Toeplitz form. Bellow this will allow us to apply 
known mathematical results concerning the asymptotic 
properties of its determinant A^v in the limit of large TV. 
Physically, this corresponds to the regime of long time r, 
i.e. to infrared asymptotics of correlation functions under 
interest. For arbitrary times r Eqs. (|38|) . (|37|) can be 
directly used for numerical evaluation of the determinant 
Ajv[/]- 

The following important point is also worth stressing. 
The Green functions of LL tunneling spectroscopy and 
FES problems require evaluation of such determinants 
A(5) at phases 5 that are not small. For example, in the 
case of relatively weak LL interaction (K close to unity) 
or small phase shift Sq for the scattering on the core hole 
in the FES problem, one needs to know the determinant 
at 6 close to 2ir, see Table Q] For strong LL interaction or 
large phase shift Sq the value of 5 can be, in principle, ar- 
bitrarily large. As was discussed in our papers^ the deter- 
minant A(5) should be then understood as analytically 
continued from the region of small d. We emphasize now 
that the present regularization (discretization of time and 
introduction of the phase factor ensuring periodicity in 
energy) implements the required analytic continuation. 
Indeed, in the original form of the determinant, Eq. Q35p. 
the information about the integer part of S/2tt was not 
explicit, which made the analytic continuation necessary. 
On the other hand, in the present regularization the inte- 
ger part of 8/2tt enters explicitly via the last phase factor 
in Eq. (|37|) . As we will demonstrate below, this allows one 
to directly compute the determinant at arbitrary large 8. 



B. Simplified analysis of asymptotics via Szego 
formula 



The long time behavior of determinants of Toeplitz 
matrices can be found using the Szego theorem and its 
extension known as Fisher-Hartwig conjecture. The con- 
dition of applicability of Szego theorem requires that f(z) 
is a sufficiently smooth function. This condition is not 
fulfilled in our case. Indeed, already at equilibrium (and 
at zero temperature) /(e) has a jump at the Fermi en- 
ergy. In non-equilibrium situations we are interested in 



/(e) will have two ("double-step distribution") or more 
such jumps. We will see, however, that the Szego for- 
mula nevertheless yields correctly the main ingredients 
of the result (dephasing rate and modified power-law ex- 
ponents). A more accurate treatment will be performed 
below in Sec. IIII CI in the framework of Fisher-Hartwig 
formula. 

Here we assume \S\ < tt. (The accurate consideration 
in Sec. IIII 01 will be performed for arbitrary 8.) We first 
consider a simple case of thermal equilibrium, when the 
determinant (|36|) can be calculated explicitly 



A T (S) 



1 



/ ttTt 



(1 



A2 t 2)5(^f) 2 V sinh7r ' rT 



(5/27T) 



(39) 



Note that the precise behavior of the functional determi- 
nant at the ultraviolet scale, t ~ A -1 , depends on the 
regularization procedure. Equation (|39p corresponds to 
a smooth cut-off e~' c '/ A in the energy spaced that is dif- 
ferent from the regularization we use in this work. We 
are interested, however, in energy scales much less than 
A. i.e. r ^> A -1 , where the determinant does not depend 
on the regularization scheme, up to an overall pref actor 
independent on the distribution function n(e). At T = 
one readily finds from Eq. (|3"9"|) 



A T (<5) ~ (rA)-W 27r)2 , r>A- x 



(40) 



To apply the Szego formula (sec Appendix [B]), we have 
to calculate the Fourier transform V(tj) of In /(e), 



V(t j ) = f A ^~ let3 In /(e). 



(41) 



where /(e) is given by Eq. ([57)1 . For the case of the 
Fermi-Dirac distribution with T = one finds 



(42) 



According to the (strong) Szego theorem [Eq. (|B4[) ]. the 
large- TV behavior of the determinant reads, in the present 
notations 



JV 

Ajy - exp <{ N AtV(0) + At ^tjVit^V^-tj 

3=0 

(43) 

The first term in the exponent TV At V(0) = tV(0) = 
— i<5rA/27r is purely imaginary and yields just a phase 
factor. Calculating the second term, we find 

N ( a \ 2 r T At 



3=0 



2tt j ./a. -1 t 



-) ln(rA) 



(44) 



7 



Here we assumed that the time is sufficiently long (com- 
pared to the ultraviolet scale), rA ^> 1, which is ex- 
actly the condition of applicability of the Szcgo theo- 
rem (N ^ 1). Exponentiation of Eq. (|4"4")) according to 
Eq. (J43]) reproduces the result (|40)) . 

It is easy to verify that the Szego formula yields the 
correct behavior of the determinant also at finite temper- 
ature T. (In fact, in the long-time regime tT 3> 1 the 
Szcgo theorem becomes rigorously applicable.) 

We turn now to the non-equilibrium situation and fo- 
cus on a double step distribution function, 



n(e) = (1 - a)n Q (e - e ) + an Q (e - ei) 



(45) 



Here eo = —all, t\ = (1 — a)U, and no(e) is the zero- 
temperature Fermi-Dirac function, no(e) = 6{— e). 
The function /(e) has now the form 

/(e) = e"^* 

e~ u , -A < e < e 

x { 1 + ( e - lS - 1) a , e <e<ei (46) 
1 , ei < e < A. 

The Fourier transform of In/, Eq. (|4"Tj) . reads 

r -i(^, e -«ot + J 8 ie -fei*) , t.^o 

^fe) = 3 (47) 
where we defined 



2tt 



In 



1 



-j<5 



1) 



(48) 



Now we apply the Szego theorem. The leading, linear- 
in-£, term in the asymptotics of In A is governed by V(0) 
given by the second line in Eq. (|47|) . Particularly impor- 
tant is the real part of V(0) that leads to the exponential 
decay of the determinat with time, A oc exp(— t/2t^). 
The corresponding decay rate is given by 



2ReV{0) = - — In 

Z7T 



1 - 4a(l - a) sin 2 - 



(49) 



To find the subleading term in the Szego formula (|4"5)l , we 
have to evaluate the sum /At^^^t^V {t k )V {—tk) . We 
get 



N 

AtJ2tkV{t k )V{-t k ) 



(50) 



k=l 



A- 1 1 



/? 2 + 2/3 /3i cosC/i) 



(/?o 2 



3 2 )lnC/T , 



T » [/" 



(51) 



In the short-time regime, Ut <C 1, we simply reproduce 
the equilibrium result. On the other hand, for long times, 
Ut ^> 1, a different behavior emerges, 



A T (6) 



~t/2t 4> 



(A/£0 



-05/27T) 



(Ut) 



(52) 



where /3 , /?i are given by Eq. (|4^)) and t$ by Eq. (|4*9")l . 

As has been already mentioned, the present problem 
goes, strictly speaking, beyond the range of applicability 
of the Szego theorem, since the function /(e) has discon- 
tinuities. This results in the (correct) \i\N behavior of 
the second term in the exponent of (|43|) . while it should 
have a constant limit as N — > oo under the conditions of 
applicability of the Szego theorem. It turns out, however, 
that the key results obtained above — dephasing rate and 
modified power-law exponents — are correct. This will be 
shown in Sec. 1111 Cl bv using recent mathematical results 
on Fisher- Hartwig conjecture which treats Toeplitz deter- 
minants of exactly the type we have encountered. The 
application of the Fishcr-Hartwig formalism will allow us 
not just to confirm the above results but also to go con- 
siderably further. First, we will calculate the asymptotics 
of the determinants exactly, including prefactors ~ A^ . 
Second, we will obtain results for an arbitrary phase 5. 
As we have already emphasized, this is important for the 
analysis of the many-body problems considered above. 
Third, we will obtain not only the leading contribution 
but also subleading terms. As we show below, various 
contributions have very transparent physical meaning in 
the problems of FES and tunneling spectroscopy, corre- 
sponding to power-law behavior at multiple Fermi edges. 



C. Accurate analysis of asymptotics via 
Fisher-Hartwig conjecture 

It is instructive to begin again by considering the zero- 
temperature equilibrium case, n(e) = 9(—e). In the case 
of a single singular point, the Fisher-Hartwig generating 
function f(z) has the form 



A) 



f(z) = e v ^\z-ZoF° ^-j Wb(*), (53) 

where g Zo ,p (z) is a "jump function", Eq. (|B6|) . Compar- 
ing Eq. (|53|) with Eq. (J3TJ) (where we set 6 ~ ire/ A and 
z = e 1 ®), we identify the parameters: zq = 1, a® = 0, 
/? = -<ty27r, and V{z) = -iS/2. Using the Fisher- 
Hartwig formula (|B8[) . we obtain the asymptotics 



A N (5) = e-^N-^ 2 G[l-^- 



G 



S 



(54) 

where G(z) is the Barnes G-function. In physical nota- 
tions N = tA/tt. It is now easy to see that the Fisher- 
Hartwig formula (|B8[) reproduces correctly the power-law 
behavior (j40|) of the determinant. 

The result (|54l) is valid for arbitrary 5. Consider an 
important case of 5 in the vicinity of 2tt, S = 2tx + 5' 
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with 1 8' | <C 1, that is relevant to the FES problem with 
small scattering phase and to tunneling spectroscopy of 
LL with weak interaction. Using G(l) = 1 and G(z) ~ z 
for small z, we get 



. " -irA(l+<5727r) I 7 A 



Z7T \ 7T 



-(1+57277) 2 



(55) 



We have introduced here /3 and /3[ satisfying |Re/3^| < 
1/2, P' + f3[ = -S'/2it. The exponents j3 3 (that may dif- 
fer from /?'■ by an integer only) satisfy /?o + /?i = — 5/2ir = 
—M — 8' /2tt. Equations (joTj) represent one possi- 

ble choice; the final result will not depend on a particu- 
lar choice in view of the summation over integers rij in 
Eq. JB9]). We obtain from Eq. (|B9|) : 



It was shown earlier— that A T (5) at 8 — > 2ir should yield, 
up to a proportionality factor, the free fermion Green 
function G (r) , sec Eq. (TTJ) . As we see from Eq. ([5~4l , the 
exact correspondence in the present ultraviolet regular- 
ization of the determinant is 



G|(r)=e- A 



— t * lim ^r(2n + S>) . (56) 



We will demonstrate below how this general formula 
works for the case of a multiple-step distribution. 

We now turn to the case of a double-step distribution 
(|4"5)) . The singular points are zj = e l7Te j/ A (j = 0,1), 
where eo = —all, e.\ = eo + U. We consider first the case 
of \S\ < 7r. It is easy to see that Eq. (|46|) is of Fisher- 
Hartwig form of generating function with two jump-type 
singularities, 

f(z) = e v ^ (£^° (±y X g Z0 ,? (z)g zlA (z), (57) 
with fa (j = 0, 1) given by Eq. gSJ and 



... . —id .Utt ( aS 

V{z) = const = — +t —(- +f)l 



(58) 



The logarithm in Eq. (|48[) and in analogous formulas for 
0'j below is understood in the sense of its main branch 
(with imaginary part between — n and w). We further 
note that under the condition U <C A we can approximate 
\zi — z \ ~ ttII/A. Applying Eq. (|B8|) . we thus get 

iS 

A T (<5) ~ exp{-— tA - ir/i - t/2tA 



A 



-2/3 /3i 



rA\- p °-^ 

x G(l + A,)G(l-/3o)G(l + /3i)G(l-/3i), 

(59) 

where l/r^ is the exponential decay (dephasing) rate, 
Eq. (09]), and /x = -E/(Re/3i + a8/2ir). This confirms 
the long-time ([/t» 1) behavior of the determinant ob- 
tained above from the Szego formula, Eq. (|52|) , and yields 
the exact value of the corresponding prefactor. 

Consider now the more general situation, when 6 is not 
small: 8 = 2ttM + 8' with some integer M and \8'\ < w. 
The exponents /3q and j3i can be now chosen to be 



Pi = — — — ln(l — a + ae 
2tt 



-iS\ 



(60) 



8' 



2 71 



-M - — - ft = -M + p' .(61) 

Z7T 



A T (8) 



id 

exp{— — rA — iTfjf — iMeoT — t/2ta,} 
2tt 



oo 

E 

n— — oo 



rA 

7T 



A \ -2(0£-M+n)(#-n) 

x e~ mUT G(l +/3' -M + n)G(l — Q + M — n) 
x G(l + # - n)G(l - # + n) , 

(62) 

where // = -C/(Re/3J +o<J727r). When M = 0, the dom- 
inant term (that is characterized by the smallest power- 
law exponent Re[(/3 — M + n) 2 + (@[ — n) 2 ]) is the one 
with n = 0, reproducing Eq. (|59[) . 

In the particularly interesting case of M = 1, when 5 
is in the vicinity of 2ir, the two leading terms are those 
with n 
find 



and n = 1. Retaining only these terms, we 



A T ((5) ~ exp{- 



tA 

7T 



j<5 

2^ 



rA — ir/i' 

K) 2 



(1 - a)e 



iaJJr 



ae i(a-l)U, [ tA 

7T 



rA 

2^; 



- T/2r } 



2^ 



(63) 



In the limit 8' -> 0, substituting Eq. jB3} into Eq. ([55]l, 
we correctly reproduce the free-electron Green function 
for the double-step distribution, g|(t) = [(1 -a)e taUr + 
ae^ 1_a ^ c/T ]G ( ^ T=0 (r), where GJ r=0 (r) is the equilibrium, 

T = value of G|(r). 

These results can be generalized to a multi-step dis- 
tribution, Fig. [T] Consider a distribution function of the 
form 



n(e) 



1 = a , 
ai , 

= a m +i 



e < eo 

e < e < ei 



(64) 



-i < e < 
< e. 



Here all Oj with j = 1, . . . , m satisfy < cij < 1 with- 
out any further restrictions. In particular, no require- 
ment of monotonicity is imposed: the distribution n(e) 
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can describe inversion of population in some regions of 
energy, with Oj+i > ctj. Using Eq. (|B9|) . we obtain for 
S = 2wM + 5' where, as before, M is the integer closest 
to 8/2n: 

iS 

A T (6) ~ exp{- — rA-irfi' - r/2r^,} 

Z7T 

n + ...+n m = — M j 
A 



n 

j<k 



J G(l + -ft) 



(65) 



Here the exponents ^ (satisfying |Re/3^ | < 1/2) are 



# = -— [ln(l - a t + a t e- ld ) - ln(l - a l+l + a l+1 e- lS )] 

(66) 

the dephasing rate reads 



z = 2Im E fa 



1 m S 
-7^ - c J-i)ln[l - 4oj(l - a j) sin2 gl 



(67) 



U jk = \ej - e k \, and p! = -Re^V P'j£j- 

The remarkable periodicity in the dependence of the 
dephasing rate for a multi-step distribution, as a function 
on the phase 5, should be emphasized. When applied to 
LL spectroscopy, this results in periodic dependence of 
I/t^ on the interaction strength^. The dephasing rate 
manifests itself in a broadening of the singularities in the 
energy space (see Sec. IIIID|) . as well as in an exponen- 
tial damping of Aharonov-Bohm oscillations in out-of- 
cquilibrium intcrfcrometry ' . 

While the above results are obtained in time represen- 
tation, the experimental measurements of Green func- 
tions describing the FES, Eq. (|4| and the LL tunneling 
spectroscopy, Eq. ([5]), are normally performed in energy 
space. It is thus important to see what the implications 
of the above findings in the energy representation are. 



D. Singularities in energy representation 

After Fourier transforming from time into energy 
space, Eq. (|65|) yields multiple power-law singularities of 

the type |e + J^j n j e j ~ p'\~ 1+ ^-' j ^ , where (3j = fi'- +rij. 
Positions of the singularities are given by linear combina- 
tions of the singular points ej of the distribution function, 
with a small overall shift //. All singularities are broad- 
ened by the dephasing rate 1/2T0. 



As an important example, consider the case M = 1, 
where S is close to 2tt. The leading singularities then 
correspond to all nj being equal to zero except for one, 
rifc = — 1. The position of each such singularity is close to 
e = e k and the exponent is J^j Pj = 1 ~ ^P'k + J2j(P'j) 2 - 
There is such a singularity for each of the singular points 
of the original distribution function (i.e. for each k = 
0, 1, ... , to) as expected. Further singularities arc much 
weaker. The next ones correspond to all rij being zero 



except for n k 
is located at e 



T and 



p 



1. Such a singularity 



ni 

efc + e/ — £pj with an exponent J2j Pj = 
3-2/3^-2/3 ; '+2/?^+X; : , (^) 2 - The next ones are generated 
by rife = —2, n; = 1, located at e = 2e k — (i and are 
characterized by exponents J2j Pj = 5 — 4/3£ + 2(j[ + 
^2j(P'j) 2 , and so on. These subleading singularities can 
be understood as resulting from inelastic processes. For 
example, the edge e = efc+e; — e p results from a creation of 
an electron near e& accompanied by creation of a particle- 
hole pair with energies near e; and e p , respectively. 

Let us present results for dominant singularities for the 
FES Green function ([¥]) in the energy space in an explicit 
form. (The LL tunneling spectroscopy Green function is 
analyzed in the same way.) The behavior of the Green 
functions (e) for energies close to singular points e k + 
// is given by (up to an additive contribution that can 
be considered as constant near e k + //) 



G*(e) 



A 



2/3; 



n 

3<l 



x (e - e k - // - «/2t P 
1 T (afe - Ofc+i) , 



A 



e < (-k 



-2^ ft' 



M (-go 1 ) 

1 ± (oft - ftfc+i) , e k + fi'<e, 



where the exponents jk are given by 7fe = ~2f3 k + 

Ej(P'j) 2 - 

Let us assume for simplicity that all distances be- 
tween consecutive singular points are of the same order, 
Ujj-i ~ U. The region of validity of the behavior (l68l) is 
then |e — e/j — < [/. Let us emphasize that the power- 
law singularity in Eq. (|68[) is smeared by the dephasing 
1/2t0. In a generic situation, when the phase 5' = 5 — 2ir 
is of order ir (i.e. not small), the dephasing rate is of 
order U. Then the smearing is strong, and the power law 
essentially does not have room to develop. On the other 
hand, when the phase <5 = 27r + <5' is close to 2ir (which 
corresponds to a weak interaction in LL or to small phase 
shift for scattering on core hole in the FES problem), the 
dephasing rate is small as 1/2t0 ~ {8') 2 U , which yields 
a parametrically broad interval for the power-law behav- 
ior, {5') 2 U < [e — ejb — fi'\ < U. Note, though, that the 
power-law exponents 7; in this situation are also small, 
7i ~ 5' ', so that the power law essentially reduces to a 
logarithmic correction. 

When applied to the problem of split FES, see Eq. (|4|), 
our result agrees with that obtained by Abanin and 
Levito\*i£ for a double-step distribution. It should be 
emphasized that our approach is different from the one 
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FIG. 1: Distribution functions with multiple edges. Two 
examples of three-step distributions are shown: (A) with 
monotonously decreasing occupation; (B) with population in- 
version. 



1 








Y 

'eq 






Y 

'eq - 














Y, ^ 


i 



used in Ref. [iH . While we work in the non-equilibrium 
bosonization framework and present the Green func- 
tion in terms of a single determinant A (5) at the phase 
S = 2it — 25o, Abanin and Levitov used the fermionic FES 
theory and obtained the result in the form of a product of 
a determinant A(— 25$) and a Green function, and then 
analyzed both terms by an approximate solution of the 
corresponding Ricmann-Hilbert problem. 

Applying these results to the LL tunneling spec- 
troscopy, sec Eq. we obtain split power-law singulari- 
ties, with modified (compared to the equilibrium regime) 
exponents and with broadening by the non-equilibrium 
dephasing rate \)2t^. This is illustrated in Fig. [5] where 
we show the behavior of the tunneling density of states 
(TDOS), i/(e) = [G<(e) - G>(e)]/2m. As discussed 
above, for a weak interaction the non-equilibrium power 
laws reduce to logarithmic corrections and are weakly 
smeared (by small I/t^). For an arbitrary strength of in- 
teraction the scale for smearing of singularities becomes 
comparable to the distance between the singular points. 
The profile of TDOS in a general situation can be ob- 
tained by numerical evaluation of the Toeplitz determi- 
nant. 

It is worth mentioning some further recent works that 
addressed the non-equilibrium LL spectroscopy. Influ- 
ence of non-equilibrium conditions on exponents was also 
found within the functional renormalization group ap- 
proach (by using approximations justified for weak inter- 
action) in Ref. 1 19c dephasing was discarded there. Qual- 
itatively similar results (modification of exponents and 
oscillatory dependence of dephasing rate on the interac- 
tion strength) were also obtained in Ref. H^. The differ- 
ence is due to the fact that the non-equilibrium setup of 
Ref. [23] (a biased quantum wire with an impurity inside 
the interacting region) is different from that of our work, 
where the non-equilibrium distribution is assumed to be 
formed by scattering outside of the interacting LL region. 

If the phase S is close to 4-7T, M = 2, the leading singu- 
larities are given by iik =n\ = — 1 , located at e = efe + (4 , 
with an exponent 0] = 2 - 2f3' k - 2(3[ + £\ (ft ■ ) 2 . The 
next singularities are produced by = —2, located at at 
e = 2e k , with an exponent £\ 0? = 4 - 4(3[ + Y.jiftf, 
etc. For the case of a doublc-step distribution and S = 4tt 




FIG. 2: Schematic results for TDOS in a LL with not-too- 
strong interaction and with multiple-step distributions. The 
distributions of electrons from both reservoirs are assumed 
to be equal and of the type shown in Fig. [TJ\ (upper pane) 
or Fig. |TJ3 (lower panel). The exponents 74 characterizing 
ZBA at multiple edges are indicated. All singularities are 
broadened by the non-equilibrium dephasing rate 1/2^. 

this agrees with the exact result for A r obtained in Ref. 0- 



IV. CONCLUSION AND 
ACKNOWLEDGMENTS 

To conclude, we have considered several many-body 
problems out of equilibrium, including the FES, the 
counting statistics, and the tunneling spectroscopy in LL 
of fermions as well as of bosons with hard-core repulsion. 
We have shown that the correlation function in all these 
problems may be expressed in terms of Fredholm deter- 
minants of counting operators. The operators are con- 
trolled by the (non-equilibrium) distribution function, as 
well as by the value of the scattering phase depending on 
the interaction strength. Our non-equilibrium bosoniza- 
tion approach allows to solve these problems and to es- 
tablish connections among them. 

We have performed an analysis of the long-time asymp- 
totics of the relevant Fredholm determinants (which are 
of Toeplitz form) . In the interesting case of double-step 
(or, more generally, multiple-step) distribution functions 
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the corresponding generating functions possess Fisher- 
Hartwig singularities induced by Fermi edges. When 
transformed from time into energy representation, the 
results reveal power-law behavior, associated with multi- 
particle processes at various discontinuities of distribu- 
tion function (edges) . The power laws thus obtained dif- 
fer from equilibrium one; in addition, the singular behav- 
ior is broadened by non-equilibrium dephasing rate. 

This work was supported by GIF, DFG Center for 
Functional Nanostructures, Einstein Minerva Center, 
US-Israel BSF, ISF, Minerva Foundation, DFG SPP 
1285, EU GEOMDISS, and Rosnauka 02.740.11.5072. 



Similar calculation for the G < component yield 



ttT I 2A 



1 



2ttT J 



exp 



2TJ 
(A6) 

The ratio between the Green functions G > and G < is 
equal to 



G<(«) 



exp- 



(A7) 



in agreement with the Fluctuation-Dissipation theorem. 



Appendix A: Hard-core bosons at equilibrium 

At equilibrium the action is quadratic and to find 
Green function one needs to calculate Gaussian func- 
tional integrals 



G%(t) = -i Po J 2^expUS[«^] 

x(l-e- iur ) + 4> L (w,q)(l + e- iur ) 
The action at equilibrium is given by 



2V2 



E 



<t>h{u,q) 



(L O R) 



S\ 



>L = 



(ott 



(f> L (-uj,-q)U a L (u),q)(j> L (w,q) 

(Al) 



+ -^i(-W,-g)II L (uj,q)4> L (uJ,q) 



Performing the Gaussian integration over bosonic fields, 
one finds the Green functions 



ITT 

~2 



2 roo 



(du))(dq) 



(A2) 



IP£ (w, q)(l — coswr) — 2i sin ujtW l (lo, q) + (L <-» R) 



The integrals over momentum q can be easily calculated 



(dq) 



2 nf(w,g) 



1 UJ 

coth — — 

27rw 2T 



(dq)q- 2 W L {u J ,q) = - 



(A3) 



Performing the standard integrals over uj one obtains 

g- 2 ? sgnr 



G b {t) = -ipo 



itTt 



sinhTrTr (l + AV 2 ) 1 /*' 
Switching to the frequency domain one obtains 



(A4) 



,1 !W \ 



2ttT i 



exp 



(-) ■ 

(A5) 



Appendix B: Mathematical background: Szego and 
Fisher-Hartwig formulas 

For completeness we present here the short summary 
of known mathematical results concerning the theory of 
Toeplitz determinants, see Ref.[58|and references therein. 
A Toeplitz matrix {fij}, ) < i, j < N— 1 is generated by a 
complex- valued function f(z) on the unit circle z = e l6 , 
where 9 is the polar angle 9 € [0, 27r]. Entries of the 
Toeplitz matrix are the Fourier coefficients 



(Bl) 



The Szego theorem is formulated for the case when 



/(*) 



is non-zero and sufficiently smooth on the 



unit circle. The function V(z) may be described by its 
Fourier harmonics 



k— — oo 



where 



v k = 



(B2) 



(B3) 



The smoothness condition requires that X)fcL-oo l^ll^fcl 2 
converges. It is further assumed the aigV(z) returns to 
its original value (rather than picking up a 27rn contri- 
bution with non-zero integer n) when z goes around the 
unit circle. According to the (strong) Szego theorem, the 
large- N asymptotic behavior of the determinant Ajv[/] 
of the corresponding matrix is 



A N [f] = exp iNVo + J2 kv kV-k 



(B4) 



fc=i 



The Fisher-Hartwig formula deals with Toeplitz ma- 
trices of a more general form, with generating function 
f(z) having m + 1 singularities (m = 0, 1, 2, . . .), 



f(z) = e v ^z^o^Y[ 



(B5) 



12 



where 




-7r < argz < 6j 



< argz < 7r . 



(B6) 



where G(z) is the Barnes G-function. A more general 
results that is valid for any values of fij and yields also 
subleading contributions has the following form: 



The singularities are located at points Zj = e j with 
j = 0, . . . , m; for definiteness, they can be assumed to be 
ordered as follows: 



- 7T = O < 



< 



< 



< IT. 



(B7) 



The strength of singularities is controlled by a set of pa- 
rameters (X,-, j3j, satisfying Re ctj > — |, /3j € C. 

Derivation of the asymptotic behavior of Toeplitz de- 
terminant with Hartwig-Fisher singularities as well as 
overview of previous literature can be found in the re- 
cent work£&. In this paper, we are interested in a par- 
ticular case of ay = and V(z) = const = Vq. In- 
deed, the function (|37[) . with distribution n(e) having 
double-stcp or multiple-step form (superposition of two 
or more zero-temperature Fermi distributions with dif- 
ferent chemical potentials) belongs exactly to this class 
of function. We thus present the results for this particu- 
lar case only, referring the reader to Ref. [H for general 
results. If all (5j are sufficiently close to each other, such 
that |Re/3j — Re/3fc| < 1 for all j, k = 0, . . . , m, the asymp- 
totic behavior of the determinant reads 

A N = e NV °N-^T=o% Y[ \zj-z k \ 20i0k 

0<j<k<m 



(B8) 



A w = e 



NV 



e n 

n + ...+n m =0 j'=0 



,N 



0<j<k<m 



in 



x [] G{l + ^)G{\-Pi) 



(B9) 



The summation in cq. (|B9[) goes over all sets of inte- 
gers no, rii, ■ ■ ■ j n m satisfying X]7Lo n i = ^- This formula 
plays a central role in our analysis of the asymptotics 
of determinants governing the Green functions of many- 
body problems in Sec. lIII Cl As we show there, the shifts 
by integers rij in Eq. (|B9[) generate contributions corre- 
sponding to multiple Fermi edges to these Green func- 
tions. 
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